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Pairing Correlations in a Generalized Hubbard Model for the Cuprates. 
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Using numerical diagonalization of a 4x4 cluster, we calculate on-site s, extended s and d^2_y2 
pairing correlation functions (PCF) in an effective generalized Hubbard model for the cuprates, with 
nearest-neighbor correlated hopping and next nearest-neighbor hopping t' . The vertex contributions 
(VC) to the PCF are significantly enhanced, relative to the t — t' — U model. The behavior of the 
PCF and their VC, and signatures of anomalous flux quantization, indicate superconductivity in 
the d-wave channel for moderate doping and in the s-wave channel for high doping and small U. 



Since the discovery of high temperature superconduc- 
tivity, much effort has been devoted to study the prop- 
erties of the Hubbard model, the t — J model, and mod- 
ifications of them. While these studies have helped to 
clarify several optical and magnetic properties of the 
cuprates HJ^], the superconducting mechanism remains 
unclear. Studies in generalized t — J models suggest a 
magnetic origin of superconductivity ||3|-|7|, but the nu- 
merical results seem to require either a superexchange J, 
or a three-site term [||J3l, which is beyond the realistic 
range for the cuprates. In addition, the constraint of no 
double occupancy in these models reduces the mobility 
of the superconducting pairs |8|. On the other hand, the 
search for signals of superconductivity in the Hubbard 
model have been negative so far ||l|,y,|l^. This fact stim- 
ulates the study of modifications of the Hubbard model 
which represent more closely the physics of the cuprates 
[p|,pd| . Recently an effective modified Hubbard model for 
the cuprates derived earlier [Q has been studied [ pT[ . 
The model includes a nearest-neighbor (NN) correlated 
hopping which depends on the occupation of the two 
sites involved and next-NN hopping t' . Within a mean- 
field approximation [ pT| , the correlated hopping has been 
found to originate pairing, the underlying mechanism be- 
ing similar to that provided by a superexchange coupling 
J P,0,Q3[ . The shape of the Fermi surface and the posi- 
tions of the van Hove singularities (vHS), modified with 
t' , influence the magnitude and the symmetry of the or- 
der parameter. The expected instability for moderate 
dopings, is d-wave superconductivity in concurrence with 
long-range antiferromagnetism near half-filling. 

In this Letter we report results on pairing correlation 
functions (PCF) and spin correlation functions for this 
effective model, obtained by numerical diagonalization 
of a square cluster containing L — 16 unit cells. We 
find evidence of strong superconducting correlations with 
d^2_y2--wsive symmetry in the doping regime of interest 
for the cuprates. Furthermore, in contrast to the case 
of the ordinary Hubbard model, we find indications of 



anomalous flux quantization (AFQ), characteristic of su- 
perconductivity jlj], in most of the explored region of 
parameters. Our numerical results support the mean- 
field picture. The size of the cluster and the inclusion 
of t' made the calculation particularly difficult. To our 
knowledge, even with t' = Q, no exact PCF have been 
so far reported in this cluster allowing doubly occupied 
sites. 

The effective model for the cuprates is [O : 
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where < ij > {< ij >) denotes NN (next-NN) positions 
of the lattice. U represents the cost in energy of con- 
structing a Zhang-Rice singlet from two singly occupied 
cells. tAA represents the hopping of a Zhang-Rice singlet 
to a singly occupied NN cell. The terms with amplitude 
tAB correspond to the destruction of a Zhang-Rice singlet 
and a nearest-neighbor cell without holes, creating two 
singly occupied cells and vice versa, i^s describes the 
movement of an isolated hole. While U lies between 3 and 
4eV, the magnitude of the correlated hopping terms is ten 
times smaller, and tAB ^ 10% larger than {tAA +tBB)/'^ 
has been estimated |l2] . However for other parameters of 
the multiband model, this ratio can be much larger, since 
tAB is linear in the Cu-0 hopping tpd, while tAA, tBB ^ 
^pd liil' ^^ mean field, for tAB > tAA, tBB, supercon- 
ductivity in the s- and d-wave channels is obtained pi| . 
Near half filling, d-wave superconductivity competes with 
long-range antiferromagnetism. If i' = 0, a SDW takes 
place at n = 1 while finite t' destroys perfect nesting, 
and for doping such that vHS lie near the Fermi level, a 
d-wave superconductor coexisting with short-range anti- 
ferromagnetic fluctuations is expected. Instead, s-wave 
superconductivity develops for small U and sufficiently 
small particle densities n. Although vHS are not well 
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defined in a small cluster, t' introduces changes in the 
distribution of the particles in fc-space and conclusions 
concerning the tendencies in the behavior of the PCF 
can be extracted. We restrict to the electron-hole sym- 
metric case tAA = tBB = 1 and large Iab ^ 2, in order 
to render more noticeably the effects of the correlated 
hopping. We also investigate t' — 0, —0.45 |l5[ . 
The PCF are: 



P^{^) = (Ai(.)A„(0)), 



(2) 



where for on-site s pairing Atg(i) — cLcL , while Aj,(i) = 

EsUsU 



^i+S^^'i.i 



h5xcI|]/%/8, with fes{S) = 1 for ex- 
tended s pairing, and fdiS) — 1 {fd{5) ~ —1 ) when 
(5 = ±(1, 0) ((5 = ±(0, 1)) for dx2-y2 pairing. We normal- 
ize Ajj(i) in such a way that |A]j(?)|0)p = 1, to facilitate 
comparison among the different PCF Mm. To compute 
the VC to the PCF 0, denoted as Pa{i), the quantity 
{{c'^c^) (cj^c^) — (c'^Ci,) (cj^c^)) is subtracted, for every term 
in Eq. (H) of the form {c\c^ Ci,c^) . For a BCS ground 
state, Pa{i) is positive and proportional to the square 
of the order parameter. The results we show for the 
correlations functions correspond to "optimum" bound- 
ary conditions (OBC), which could be periodic (PBC), 
antiperiodic (ABC) or mixed (MBC), i.e. periodic in 
one direction and antiperiodic in the other, according to 
those which lead to the minimum ground-state energy. 
The computation has been made possible by exploiting 
all symmetry operations of the space group of the square 
lattice [Q plus time reversal (256 operations in the clus- 
ter). Half of these operations are lost for MBC, and the 
reported PCF are averages over equivalent distances in 
the periodic system. 

To give an idea of the expected magnitude of the PCF 
and in order to establish a criterion to interpret our re- 
sults, we analyze the behavior of the PCF and the VC for 
the usual attractive Hubbard model with a quite large at- 
traction U = —5, in which case superconductivity is well 
supported by several calculations ||l^. These quantities 
are displayed in Fig. 1 for A^ = 10 particles and distances 
larger than one lattice site ^y] . As in the case of previous 
Monte Carlo results |1^, Pair) shows oscillations with 
distance r, while Pair) exhibits a smoother behavior. It 
is clear that Pos{r) dominates over the other PCF, which 
is in agreement with the s-wave, predominantly on-site 
character of the superconductivity in the model |lSf| . In 
the light of these results, we establish the following cri- 
terion to extract information from our numerical data: 
we conclude that superconducting correlations in the a- 
channel are present in the model when both quantities, 
Pa and Pa, are enhanced at large distances relative to 
the non-interacting case. 

In Fig. 2, we show the effect of tAB and t' for U — 0. 
The PCF Pa{r) (not shown) display the same qualita- 
tive behavior as those in Fig. 1 (a). We conclude that 



for these parameters the model has strong signals of s- 
wave superconductivity in both on-site and NN channels. 
This agrees with the mean-field calculations |1^. For 
N = 12 particles, the values of Pair) (not shown) are 
reduced in ^ 0.01, but the qualitative behavior remains 
the same. For t' = 0, superconducting correlations are 
also strongest in the s-channel. For N = 12, Pes{V8) is 
approximately one half of the corresponding value for 
A^ = 10. For both densities, a negative t' enhances 
Pes,os{r) relative to the case with t' — 0. According 
to Ref. flit], when U overcomes a certain value, s-wave 
superconductivity is replaced by a SDW when t' — 
and by rf-wave superconductivity for finite t' . Keeping 
t' — —0.45, tAB = 2, and increasing U, we find a de- 
crease in the s-wave PCF and an increase in the d-wave 
ones. The latter dominate already for U = 4 and A^ = 10, 
with Pd{r) ~ 0.015 and values significantly larger than 
those for tAB — 1- 

For U = 10, with tAB — 2, the PCF are much larger in 
the d-wave channel. The behavior of Pd{r) and Pd{r) 
for different densities (A^/L=0.625, 0.75 and 0.875) is 
shown in Fig. 3, for t' = 0, —0.45. To simplify the fig- 
ure, we do not show the values of -Prf(r) for tAB = 1 
and for the non-interacting case. For A^ = 10, 12, with 
tAB = 2 and t' = -0.45, the values Pd{r) ~ 0.02, 0.03 
(Fig. 3(d)) at distances \/2 < r < ^/8 are roughly 
half of the values of Pqs (r) for the Hubbard model with 
strong on-site attraction (Fig. 1), and very similar to 
those of Pd{r) for a short-range resonance- valence-bond 
wave function which by construction has superconduct- 
ing ODLRO |H,f6| . These results are strong indications of 
d-wave superconductivity. We should also note that the 
superconducting d-wave pairs in the model, have an in- 
ternal structure which extends beyond NN and with only 
a partial overlap with A]j(i). Thus, our d-wave PCF are 
reduced with respect to the optimum normalized PCF by 
the square of this overlap [g[ . 

For A^ — 10, 12, the effect of a negative t' is to en- 
hance the VC Pd{r). Instead, for A^ = 14, both Pd{r) 
and Pd{r) are large for the case with t' = 0, while they 
are very small for t' = —0.45. Note that in all the cases 
with sizable pairing correlations, the values of Pdir) cor- 
responding to tAB = 2 are significantly larger than those 
corresponding to tAB = 1, with the same values of t' 
and U. In addition, in these cases, the non-interacting 
Pd{r) lie bellow the displayed ones for tAB = 2 in Fig. 
3 (a). The remarkable large values of Pd{r) observed 
in Fig. 3(c) for the case with A^ = 14 particles could 
be somewhat exaggerated due to particular finite-size 
effects 21|. In fact, when tAB = 1, Pd{r) in Fig. 3 



(c) is large while the values of Pd{r) are smaller than 
those of the non-interacting case. The mean-field treat- 
ment mW predicts a maximum of the superconducting 
gap with d-wave symmetry at half- filling for i' = 0, 
when long-range antiferromagnetism is not taken into ac- 
count. The concurrence between superconductivity and 
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long-range antiferromagnetism near half-filling manifests 
itself in this cluster when different BC are used. For 
N — 14, t' = 0,tAB = 2,U = 10, spin-spin correlations 
(not shown) are much stronger for PBC than for MBC, 
while in the first (latter) case Pd{r) is weaker (stronger) 
than in the non-interacting case. In any case, as expected 
P,RO|, the maximum of the PCF with doping shifts to 
higher doping as t' increases. 

In contrast to the cases without correlated hopping, 
we find signs of AFQ in most of the explored parameter 
space. AFQ consists of a periodicity of half a flux quan- 
tum in the ground-state energy E{^) as a function of a 
flux $ threading the system in a toroidal geometry, and 
it is a necessary but not sufficient condition for super- 
conductivity [IJ]. In finite systems, a tendency to AFQ 
is indicated by a crossing of energy levels with different 
total wave vector as $ is varied, and the presence of two 
relative minima in E{^) with a difference of <& in tt (usu- 
ally at <I> = and <I>=7r) in the interval [0, 27r). In Fig. 4 
we show the dependence on $ of the lowest energy levels 
of the system for several values of the parameters in which 
the level crossing occurs. Fig. 4 (a) and (b) correspond 
to dominant s-wave PCF. In the other cases shown, the 
d-wave PCF are the largest ones. For t' — Q,N — 10 we 
found a very similar behavior of -E($) to that observed 
in Fig. 4. Due to the fact that the introduction of a flux 
breaks the space group symmetry and increases consid- 
erably the size of the irreducible subspaces, we have not 
constructed E{^) curves for N > 10. 

We have also calculated charge and spin correla- 
tion functions. The spin structure factor S'(q) = 
J2^j{S■S^e''i^^^~^^^)/L^ for iV = 12 is shown in Fig. 
5. In all cases, the increase of Iab tends to restore a 
peak at (tt, tt), which is rather broad, indicating the pres- 
ence of short-range antiferromagnetism, similar to that 
found for an RVB state with superconducting ODLRO 
0. For the case N = U,t' = 0,tAB = 2, f/ = 10 5(q) 
exhibits a broad structure at (tt, tt) for MBC, while large 
Pd{r) is obtained, as discussed above. Instead, for PBC a 
much narrower peak, suggestive of longer-range antifer- 
romagnetic correlations is observed in S{q) while Pd{f) 
are weaker than those of the non-interactive case. In one 
dimension (ID), there is analytical and numerical evi- 
dence that for small U and large tAB^ the ground state 
at half filling consists of singlet dimers, and singlet PCF 
dominate when the system is doped |l3] . The natural ex- 
tension of this scenario to 2D, seems to be a short range 
RVB-like state at half filling, which turns into a singlet 
superconductor as the dimers acquire mobility with dop- 
ing. While we expect long-range antiferromagnetism in 
the half-filled case, our results for PCF and 5(q) are con- 
sistent with this scenario as the system is doped. 

In summary, we have shown that correlated hopping, 
which arises naturally in a low-energy reduction of multi- 
band models for the cuprates, leads to pairing correla- 
tions of magnitude similar to that observed in the same 



cluster for models for which superconductivity is well es- 
tablished. For low to moderate doping, the favored sym- 
metry is da;2„y2 for finite U and t' . Part of the numerical 
work was done at the Max-Planck Institut PKS. L. A 
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For i' = 0, iV = 14, the OBC are MBC. As a conse- 
quence the ground state is highly degenerate in the non- 
interacting limit, and this can affect the behavoir of the 
PCF. 




FIG. 1. (a) Pairing correlation functions and (b) vertex 
contributions to them as functions of distance for [/ = (open 
symbols) and U — —5 (solid symbols), with tAB = 1, t' — 0, 
A'' = 10 and PBC. Triangles, circles and squares correspond 
to on-site s, extended s and d^2_y2 PCF respectively. Open 
symbols coincide in (b). 



FIG. 2. Vertex contribution to the PCF for tAB = 2, 
t' = -0.45t, U = 0, N = 10 and PBC. The values without 
correlated hopping {tAB ~ 1) are also shown for comparison. 
The meaning of the different symbols is the same as in Fig. 
1. 
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FIG. 3. Pd{r) and Pd{r) for U = 10 and tAB = 1 (open 
symbols) and tAB ~ 2 (solid symbols). Circles, squares and 
triangles correspond to A = 10, 12, 14 respectively. 
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FIG. 4. Lowest relevant energies as a function of fiux for 
i' — — 0.45t, A'' = 10 and several values of U and Iab- The 
value E{0) is subtracted. Only half of the interval [0, 2n) is 
shown because E{^) = £( — $) = i5(27r — $). 



FIG. 5. Spin structure factor as a function of wave vector 
for t' = -0.45t, U ^10, N = 12. Open (sohd) circles denote 

tAB ~ l(iAS = 2). 



